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A scaling analysis is made of the order parameter describing monopole condensation at the decon- 
fining transition of Nf — 2 QCD around the chiral point. In accordance with scaling properties of 
the specihc heat, studied in a previous paper, scaling is consistent with a hrst order transition. The 
status of dual superconductivity of the vacuum as a mechanism of color confinement is reviewed. 


I. INTRODUCTION 

This paper is the fourth of a series Hii ^ on a re¬ 
search program aiming at developing and testing on the 
lattice the working hypothesis that confinement in QGD 
is produced by dual superconductivity of the vacuum, i.e. 
by condensation of magnetic charges. 

Such a mechanism implies that the deconfining tran¬ 
sition is an order-disorder phase transition, i.e. a true 
phase transition, described by an order parameter. 

Another basic underlying idea is duality ii: the 
magnetic excitations which are expected to condense are 
non local in terms of the helds. They should become local 
in the dual description, and the strong coupling regime 
of the confined phase should be mapped by duality in the 
weak coupling regime for dual helds. 

Magnetic charges exist in gauge theories [E 0 0 and 
are nicely described in a formulation of the theory based 
on parallel transport, as the Wilson formulation is 0- 

The prototype theory is the lattice 17(1) gauge theory 
in 4d. In the Wilson’s formulation this theory shows a 
deconhning transition at /3 = 2/g^ ~ 1.01, from a strong 
coupling phase in which static charges are conhned, to 
a Coulomb phase. Dual magnetic excitations can be de- 
hned and with them an order parameter which is the vev 
(/i) of a magnetically charged operator pL. This vev is dif¬ 
ferent from zero in the confined phase, and strictly zero 
in the deconfined phase. Thi s op erator is defined by its 
correlation functions [l(i fill Fl^ , exactly like the kink’s 
correlators in the Ising model of Ref. |E ■ It can however 
be explicitly defined as a non local operator in terms of 
the gauge helds [Hill in the same way as the kink cre¬ 
ation operator can be explicitly written in terms of the 
spin variables 0 

The operator /i is a gauge invariant Dirac like 00 
0 magnetic charged operator, non local but obeying 
cluster property, and a rigorous proof exists that it is the 
correct order parameter, both for the Villain 0 and 
for the Wilson ^1 actions. In the presence of electric 


^ In the following we shall refer to the previous papers as I, II, III. 


charges, however, there are problems for the existence of 
the duality transformation . We will come back below 
in Sect. ^ 

The translation of the above formalism into non 
abelian gauge theories is non trivial. The original pro¬ 
posal to dehne and expose monopoles is based on a 
procedure known as abelian projection. A gauge is hxed, 
e.g. by diagonalizing any observable in the adjoint rep¬ 
resentation, a residual gauge invariance is left 

and the monopoles are dehned as monopoles of each of 
these 17(I)s, as for the abelian case. They show up as 
singularities of the gauge transformation which brings to 
the abelian projected gauge, in the sites in which two 
eigenvalues of the chosen observable coincide. 

In this way the number and the location of monopoles 
in a given field configuration depends on the abelian pro¬ 
jection, i.e. on the choice of the gauge, and there is a 
functional inhnity of ways to do it. 

Therefore if one pretends that the monopoles detected 
on the lattice configurations are the relevant degrees 
of freedom for conhnement (monopole dominance |lf3l|). 
then one has to select the abelian projection which de¬ 
fines them. The usual choice is the so called maximum 
abelian projection, by which a gauge is selected maximiz¬ 
ing the role of the residual abelian degrees of freedom. In 
some literature this choice is called the abelian projection 
tout-court. According to this approach the monopoles in 
this projection are the effective degrees of freedom for 
confinement. 

An alternative attitude is to look at the change of sym¬ 
metries at the transition, postponing the identihcation of 
the relevant excitations, which is a far more difficult prob¬ 
lem. This amounts to define and test an order parameter 
(/i) for monopole condensate. 

The lattice version of the operator p, has been con¬ 
structed in I, II. p{x, t) creates a monopole in a given 
abelian projection at the site (x, t). p is color gauge in¬ 
variant, magnetic 17(1) gauge invariant 0|, but of course 
is abelian projection dependent. 

In III it was shown numerically that the order param¬ 
eter {p) is abelian pro jection independent. Arguments 
were given in Ref. justifying theoretically this ob¬ 
servation. Below we will present a complete argument 
showing that the statement {p) yf 0 or {p) = 0 is abelian 
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projection independent. Any /r creates a monopole in all 
abelian projections. The recent observation of the exis¬ 
tence of abelian Abrikosov flux tubes even without fixing 
the gauge fully supports this view pH. 

Creating a monopole is a better defined operation than 
detecting it; the same is true also for vortices [ 2 ^ . 
Their creation is gauge independent; their detection de¬ 
pends on the choice of the gauge. 

In summary {^) proves to be a good order parameter 
for the quenched theory: it is gauge invariant, abelian 
projection independent, strictly zero in the deconfined 
phase in the infinite volume limit [HQEi , non zero in 
the confined phase and its scaling with volume is consis¬ 
tent with the critical indexes of the transition. 

On this basis (/i) can be tentatively used as an order 
parameter for confinement also in full QCD. This would 
fill a gap: indeed in full QCD no symmetry exists, the 
chiral phase transition is defined at vriq = 0 , but an order 
parameter for confinement is missing. 

In Ref. it proved numerically that, for 

Nf = 2, (n) 7 ^ 0 below the “transition line” and {fj) = 0 
above it. In this paper we discuss a finite size scaling 
analysis of the order parameter for Nf = 2 QCD. We 
find that it is consistent with a first order phase tran¬ 
sition. The same indication comes from the analysis of 
the specific heat and of the chiral order parameter at the 
chiral transition 

In Sect. HDwe review the status of the order parameter 
(/r): its independence on the abelian projection, the diffi¬ 
culties related to the coexistence in the theory of electric 
and magnetic charges. 

In Sect. rnriwe present the numerical results for Nf = 2 
QCD and a scling analysis of the order of the deconfine¬ 
ment phase transition. 

Sect. 1^ contains the conclusions and perspectives. 


II. THE (DIS)ORDER PARAMETER {/i) 


with 


a Nc — a 



The gauge transformation U identifies the represen¬ 
tation in which (jf‘{x) is diagonal, i.e. the abelian 
projection. Of course depends on I/, i.e. on the 
abelian projection. 

In the abelian projected gauge (j)°' = which is x 

independent, the longitudinal part of E in Eq.^ cancels 
in the convolution with bj_, and only the diagonal part of 
E contributes to the trace on color indexes. If we define 
the matrices T“ (a = 1,..., A^c — 1) as 

= diag( p,..., 0 ,!,-^l, 0 ,..., 0 ) 

Nc 

these matrices are a complete basis for traceless diagonal 
matrices. Therefore 


7? _ \ ' T?a rpa 

diag — / . ^J_ diag^ 

and since 

^^{4>lagT^} = 

in the abelian projected gauge 


t) = exp 


i 



a 

_L diag 


{y,t)b±{x 



is the conjugate momentum to A“ , so that in the 
Schrodinger representation 

/r“(f,t) |A^(f,t)) = \A°^{x,t) + b±{x - y)) 


The operator jE {a = 1,..., W — 1) which creates a 
monopole of the species a in a given abelian projection 
is defined in the continuum as l28l 




exp 


i /dy Tr{(/)“(y,t)£’(y,t)}&_L(f - y) 


• ( 1 ) 


Here 


- bx{x — y) is the vector potential produced by a 
Dirac U{1) monopole sitting at x in the transverse 
gauge: 


V • &_L = 0 ; V A 6 j_ (r) 


27r r . 

- 7 -|- Dirac String; 

g r-^ 


- (/)“ transforms in the adjoint representation and has 
the form 


creates a Dirac monopole in the residual C/(l) gauge 
field in the direction of T“. 

On the other hand it can be proved dim that the 
’t Hooft tensor corresponding to a scalar field (p in the 
adjoint representation 

= Tr (3) 

reduces to an abelian form if and only if </> has the form 
(j)°‘ of Eq. [21 In that case indeed Eq. 0 becomes 

= Tr - d^irA^) - '-r 

the quadratic terms in in the definition of E^^, 

cancel and 





= 0 , 


rix)=UHx)rd^agU{x) 


( 2 ) 
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i.e. the Bianchi identities are obeyed except for sin¬ 
gularities. In the abelian projected gauge = 0, 

{A^)d^ag = and 

FY = d^Al-d,Al. 

For the simple case of SU{2) a only assumes the value 

0 = 1 , 

(l^diag = ^ (J = U\x)^a3Uix) . 

For SU (2) any hermitian field 4>{x) in the adjoint repre¬ 
sentation can be put in the form (j){x) = {x)^a 3 U{x); 

meaning that an abelian projection is obtained by diag¬ 
onalizing (/>(a:). Monopoles will appear when the gauge 
transformation U{x) is singular, and this happens when 
|(/)(a;)J _= 0 , i.e. when two eigenvalues of (j){x) coin¬ 

cide [T^ . 

For generic SU(N) any hermitian field (j){x) in the ad¬ 
joint representation can be written as 

with U{x) the gauge transformation which diagonalizes 
4>{x) and orders the eigenvalues (for instance in a decreas¬ 
ing order for their imaginary parts). U{x) will be singular 
at points where any of the coefficients Ca(x) is zero. In 
these points two eigenvalues of (j) coincide. Indeed from 
the definition of in the abelian projected gauge 

and (j)ii{x) - 0(i+i)(i+i)(a;) = Ci{x). The ma¬ 
trices (jidiag play in SU{N) the analogous role as 0 - 3/2 in 
SU{2). 

One can define magnetic currents 
j^ix) = . 

Those currents are zero everywhere, because of Bianchi 
identities, except at singular points where monopoles are 
located. In any case 

dgfg, = 0 

defining a magnetic 17(1) symmetry, and conservation 
of magnetic charge. If this symmetry is realized a la 
Wigner, i.e. if the vacuum has definite magnetic charge, 
the system is normal. The vacuum is a magnetic super¬ 
conductor if it does not have a definite magnetic charge, 
but is a superposition of states with different magnetic 
charges. (^“) 7 ^ 0 can signal the second possibility. In 
the normal state (/i“) = 0 

The lattice version of has been constructed 
in I and II, and is equally well defined in the presence 
of dynamical quarks. It was numerically shown in III 
that being (/r“) = 0 or (^“) 7 ^ 0 is a property indepen¬ 
dent of the choice of the abelian projection, and is shared 
by the version in which is put diagonal configuration 
by configuration, which is kind of an average on infinitely 
many abelian projections. The same thing results from 


comparison of measurements of (/x“) in different abelian 
projections with measurement made by use of Scrodinger 
functional as in Ref. |^. More recently flux tubes 
were observed in SU{2) gauge configurations in the 0-3 
direction without fixing the abelian projection pH| . 

In fact, if the density of monopoles is finite, i.e. if in 
whatever abelian projection there exists a finite number 
of them per unit volume, then the gauge transformation 
V(x) which connects two abelian projections is continues 
everywhere except at a finite number of points, where 
monopoles are located. 

Creating a monopole with /i“ (x) will add a singularity 
say in one of the abelian projections: the probability 
that X coincides with the position of a monopole already 
present will be zero, and there will be a neighbourhood 
of X in which V is regular. Adding a singularity in one of 
the abelian projections amounts to add a new singularity 
also in the others. creates a monopole in all abelian 
projections. Therefore (/x“) 7 ^ 0 means Higgs breaking of 
magnetic charge of type “a” in all of them. 

Dual superconductivity is an absolute property, inde¬ 
pendent of the abelian projection. 

The density of monopoles is known to be finite lU ■ A 
direct test of that can be made by looking at the distri¬ 
bution of the difference between eigenvalues of any her¬ 
mitian operator in the adjoint representation. Finding, 
in any site of the lattice, a zero difference would mean 
that a monopole is sitting on the site. The lattice can 
be made finer and finer by increasing /3. Fig. ^ shows 
that in no site the eigenvalues coincide, and hence that 
the density of monopoles is finite. The same happens at 
higher values of [3. 

A last point concerns the problem raised in Ref. 0 
about the difficulty to obey Dirac condition for mag¬ 
netic charges in the presence of electric charges. This 
is a general problem: in C/(l) this difficulty shows up 
as the non existence of the duality transformation, and 
this produces problems at short distances: an operator 
like /r should be modified at distances 0{a) by terms 
0{a^). The authors of Ref. 0 suggest a modified for¬ 
mulation in which a coulombic magnetic field produced 
by a monopole is replaced by a superposition of radial 
flux tubes with quantized flux |32|. This is a big revi¬ 
sion in the description of a monopole, but is in any case 
unmanageable from the numerical point of view |33| . In 
any case our operator is defined up to term 0{a^) as 
explained in I Sect. II, so that looking for 0{a^) correc¬ 
tions would be meaningless anyhow. Contrary to U(l) 
0{a'^) corrections in QCD should become unimportant 
in the continuum limit because of asymptotic freedom. 
The problem is however open and deserves further inves¬ 
tigation. 


III. A/ = 2 QCD 

The study of the QCD finite temperature transition 
for 2 flavors is of fundamental importance. The phase 


4 


transition is well understood at high masses (m > 2.5 
GeV), where the quarks decouple; the transition is first 
order and the Polyakov line is a good order parameter. 
At m ~ 0 a chiral transition exists, where chiral sym¬ 
metry is restored, and the chiral condensate (ipip) is an 
order parameter. At some temperature also the axial 
C/a(1 ) is expected to be restored: indeed the topological 
susceptibility drops to zero around Tc [s^. In principle 
at m ~ 0 there are 3 transitions (chiral, axial Ua(1), 
deconfinement): it is not clear if they coincide. An ef¬ 
fective description of the chiral transition can be given 
in terms of an effective free energy Assuming that 
the scalar and pseudoscalar modes are the relevant crit¬ 
ical degrees of freedom, there is no infrared stable fixed 
points for Nf > 3 and the transition is expected to be 



FIG. 1: Normalized distribution of the difference between 
the phases of consecutive eigenvalues (in units of tt) of the 
plaquette operator (upper figure) and of the Polyakov loop 
(lower figure). The values 0 or 2 correspond to two coinci¬ 
dent eigenvalues and thus to the presence of monopoles. The 
distribution has been measured at different values of /3 and 
approximately equal physical volumes, in order to check the 
continuum limit. 


first order. For Nf = 2 the transition is first order if the 
anomaly is negligible (rrijj' w 0 ) at Tc; it can be second 
order with symmetry 0(4) if the anomaly survives the 
chiral transition. In the first case the transition surface 
around m = 0 is first order. If the chiral transition is 
instead second order the surface is a crossover, and a tri- 
critical point is expected in the fJ- — T plane (see e.g. [s^), 
detectable by heavy ion experiments. 

The first scenario is compatible with dual supercon¬ 
ductivity of the vacuum being the mechanism for confine¬ 
ment both in pure gauge and in presence of dynamical 
quarks, independently of their mass: the deconfinement 
transition is then always an order-disorder phase transi¬ 
tion, not a crossover. If the second is the case instead, the 
deconfining transition cannot be order-disorder, there is 
no order parameter for confinement, and a state of a free 
quark can continuously be transfered below the “decon¬ 
fining temperature”: confined and deconfined then lose 
a definite meaning. 

The question can be answered by a numerical study 
of the chiral phase transition and has been the subject 
of extensive, even if not conclusive, literature [S^, IH, 
US ES ES ES- 113 have presented a 

large scale finite size scaling analysis of the 2 flavor chi¬ 
ral transition by looking at the specific heat, the chiral 
susceptibilty and the equation of state: we have found 
clear disagreement with a second order 0(4) behaviour 
and solid indications of the possible presence of a first 
order phase transition. On the other hand the fact that 
dual superconductivity is at work also in the presence of 
two dynamical flavours has been already sho wn, for fi¬ 
nite quark masses ~ 0.505), in Ref. [^, using 

(/j.) as a disorder parameter, and also in the Schrodinger 
functional approach in Ref. [i^ . 

In this Section we present a finite size scaling analysis 
of (^) around the chiral phase transition: if (/r) is the cor¬ 
rect (dis)order parameter for confinement, this analysis 
will give the correct critical indexes. 

The parameter (/i) is defined on the lattice as (see I 
and II) 

M = |. 

z = j (VU) , 

Z = j [VU) . (4) 

Z is obtained from Z by changing the action in the time 
slice Xq, S ^ S = S + AS. In the Abelian projected 
gauge the plaquettes 

Xlio{x,xo) = 

= Ui{x,xo)Uo{x -I- i,xo)Uj{x,xo + b)UQ{x,xo) (5) 
are changed by substituting 

Ui{x,xo) Ui{x,Xo) = C/i(f,a;o)e*^^’L(““^) 


( 6 ) 
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FIG. 2: Finite size scaling of p according to the critical indexes 
of first order transition. 

where T is the diagonal gauge group generator corre¬ 
sponding to the monopole species chosen. The numer¬ 
ical determination of (/i) is very difficult, since (/x) is 
expressed as the ratio of two partition functions. Instead 
of (/i) we measure the quantity 

(7) 

It follows from Eq. 0 that 

p={S)s-{S)s, ( 8 ) 

the subscript meaning the action by which the average is 
performed. In terms of p 

{p) = exp p(/3')d/3'^ . (9) 

A drop of (p) at the phase transition corresponds to a 
strong negative peak of p. A negative value of p diverging 
in the thermodynamical limit in the deconfined phase 
means {p) being exactly zero on that side (see I, II, III). 

We have made simulations with two degenerate flavors 
of Kogut-Susskind quarks, using the standard gauge and 
fermion actions. Configuration updating was performed 
using the standard Hybrid R algorithm ^5. The lat¬ 
tice temporal size was fixed at Nt = 4. Different spatial 
sizes {L — 12,16, 20, 24, 32) and values of the quark mass 
were used. For a more detailed account on simulation 
parameters we refer to p^. 

We can assume the following general scaling form for 
{p) around the phase transition: 

(p) = ,mLy’') 

where r is the reduced temperature and m the quark 
mass. Sending L to infinity keeping ^/L or 
fixed, must be a convergent limit. Analyticity argu¬ 
ments [ 41 L [ 4 ^ suggest that in the infinite volume limit 
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FIG. 3: Finite size scaling of p according to the 0(4) critical 
indexes. 

the mass dependence in the scaling function factorizes, 
so that P = ^ lii(/i) does not depend on the mass. In 
fact, the dependence on mLy^ must then cancel the de¬ 
pendence on the factor in front. We then obtain the 
following scaling law: 

p = L1 /‘'</.(tL1/") . (10) 

the same as in the quenched case. 

In Figure 2 we show the quality of scaling assuming 
V = 1/3, i.e. a first order phase transition: a good agree¬ 
ment is clearly visible. The deviations from scaling in 
the deconfined region are expected and are related 
to the disorder parameter {p) being exactly zero on that 
side. 

For comparison we show in Figure 3 the quality of scal¬ 
ing assuming the 0(4) critical index v = 0.75. 

The 0(4) universality class is clearly excluded, while 
there is good agreement with a first order phase transi¬ 
tion. This confirms results obtained through an analysis 
of the specific heat and of the chiral susceptibility . 


IV. CONCLUSIONS 

QCD with Nf = 2 is a specially interesting system to 
investigate confinement. The order of the chiral phase 
transition can indeed decide if the deconfining transition 
is a real transition, as required by a mechanism of confine¬ 
ment like dual superconductivity, or a trivial cross-over. 

More generally the very possibility of defining confined 
versus deconfined relies on the existence of an order pa¬ 
rameter. This parameter exists for the quenched theory, 
where the transition is first order, but is not defined for 
full QCD. For full QCD only a chiral order parameter 
('0'(/j) is defined, and only at low enough quark masses. 

We have shown that an order parameter {p) can be 
defined, which, independently of the presence of quarks 
and of their masses, characterizes confinement and de- 
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confinement, (fi) 7 ^ 0 at T < T^, {fj) strictly zero in the 
limit y —> 00 for T > T^. 

A finite size scaling analysis gives the correct 
(pseudo)critical indexes both for the quenched and for 
the Nf = 2 case. For Nf = 2 the deconfining transi¬ 
tion as seen by looking at (/i) is first order. This is the 
main result of the present article. The same indication 
comes from the analysis of the specific heat , which is 
independent of any prejudice on the mechanism of con¬ 
finement. 

Further study is on the way with improved actions and 
algorithm to settle this issue, as well as to study the 


behavior of to,;' at the chiral transition point. 

A consistent picture starts to emerge for confinement. 
A theoretical problem in the definition of (fi) in the pres¬ 
ence of electric charges exists, which deserves further 
study. 
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